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The problem of energy harvesting from flutter instabilities in flexible slender structures in axial 
flows is considered. In a recent study, we used a reduced order theoretical model of such a system to 
demonstrate the feasibility for harvesting energy from these structures. Following this preliminary 
study, we now consider a continuous fluid-structure system. Energy harvesting is modelled as 
strain-based damping and the slender structure under investigation lies in a moderate fluid loading 
range, for which the flexible structure may be destabilised by damping. The key goal of this work 
is to analyse the effect of damping distribution and intensity on the amount of energy harvested 
by the system. The numerical results indeed suggest that non-uniform damping distributions may 
significantly improve the power harvesting capacity of the system. For low damping levels, clustered 
dampers at the position of peak curvature are shown to be optimal. Conversely for higher damping, 
harvesters distributed over the whole structure are more effective. 



Increasing energy demands motivate the interest in energy harvesting concepts, where the idea is to harness the 
energy of naturally occurring phenomena. At the scale of kilowatts, concepts include energy harvesting from tidal 
currents pQ and ocean waves [2]. At the lower end of the power spectrum, concepts based on photo/thermovoltaics 
and magneto/piezoelectrics show the scope for powering sensors and mobile electronic devices [3]; these include energy 
scavenging from ambient vibrations in structures such as buildings and bridges and oscillatory motion of wheels in 
automobiles or turbines in engines [4, 5 . Energy harvesting from fluid-structure interactions (FSI) includes concepts 
such as vortex induced vibrations ( VIV) of bluff bodies in cross-flow [6j [7] , resonant vibrations induced in aerofoils 
mounted on elastic supports [8] [9], flutter of flexible plates [T0ifT2] , and combinations thereof, as for the coupled 
VIV- flutter energy harvester examined by [13]. In this work we focus on harvesting energy from flutter instabilities 
of slender structures in axial flow. 

The classical description of flutter instabilities is self-sustained oscillations that arise due to the unstable coupling of 
fluid dynamic pressure and structural bending modes, where for undamped structures the critical speed at flutter onset 
depends on fluid as well as structural properties [14]. Flutter instabilities are observed in a diversity of configurations, 
which may be broadly classified as internal or external flows [TJJ [15] . Internal flow instabilities are observed in flexible 
pipes and channels and are invariably motivated by biological phenomena such as flow induced oscillations in airways 
and veins p~6j[T7]. The pipe-conveying fluid is a canonical problem that yields deep insights into FSI; in particular 
[T8] examined the relationship between local and global instabilities and showed that locally stable configurations may 
become unstable due to wave reflections at finite boundaries; this destabilising mechanism was recently predicted for 
compliant channels as well [19 . External flow based instabilities include flapping flags [20-22] and panels [23] in a 
steady flow. For all aspect ratios, the plate can become unstable due to fluid loading, defined as the ratio of fluid and 
structure inertia, and may be represented as a nondimensional length [24] or mass ratio [25]. [26] and more recently 
[27] examined the occurrence of instabilities in slender structures. They used experimental and theoretical techniques 
to examine the role of inviscid and viscous drag on static and dynamic instabilities that arise in flexible cylinders in 
axial flow. For a recent review on the flutter dynamics of flexible bodies in external flow, the reader is referred to [28 
and references therein. 

In this work, we analyse the scope for harvesting energy in slender elastic cantilevered structures that flutter in 
a steady flow. From the point of view of the fluid-solid system, energy harvesters are essentially an energy sink, 
therefore for the theoretical approach adopted here they are modelled as internal damping in the structure. As a first 
step, we examined the feasibility of this concept using a nonlinear model of a reduced-order system, consisting of a 
slender cylinder pair connected by discrete springs and dampers [29 . It was shown that the optimal configuration for 



I. 



INTRODUCTION 



*Electronic address: kiran.singh@cantab.net Present address: OCCAM, The Mathematical Institute, 24-29 St Giles, OX13LB, Oxford 
t Electronic address: sebastien.michelin@ladhyx.polytechnique.fr 
* Electronic address: delangre@ladhyx.polytechnique.fr 



FIG. 1: Slender cantilevered structure placed in a steady axial flow of velocity Uoo^x and fixed in O. The instantaneous 
deformation of the inextensible beam is measured by the tangent angle 9(s,t) with respect to e x of the local tangent r. 

this two- degree-of- freedom system is one with energy harvesters positioned away from the instability source; such a 
configuration maintains self-sustained flapping in the presence of structural damping. 

In this work we generalise this approach for a continuous system, and seek to maximise energy harvesting through 
carefully-tailored distributions of structural damping. Conventionally, structural instabilities are stabilised by damping 
P~6j [30] . However, damping may be destabilising under moderate to heavy fluid loading conditions, as observed in 
infinite plates [31], flags [11] and fluid- conveying pipes [32]. [29] show that energy harvesting requires the presence 
of two traditionally competitive elements - flutter oscillations and damping. A configuration for which flutter is 
destabilised by damping is especially interesting from an energy harvesting perspective. This idea has been exploited 
for piezoelectric based energy harvesting from flapping flags: [TT] used linear theory to show the gain in conversion 
efficiency for flags with high fluid loading destabilised by piezoelectric based damping. In this work, we explore this 
idea using a general model for structural damping (equivalently energy harvesting) and a nonlinear model of the fluid 
structure interaction of a slender structure in a mean flow. 

It is worth noting that existing insights on damping are generally based on the assumption of a constant distribution 
of damping in the structure [see for example 10 , whilst little is known about how nonuniform damping distributions 
affect the dynamics. In this work we seek clarity on the role of damping distribution on the flutter response of slender 
structures. The specific motivation is to identify physical mechanisms that maximise this dissipated (i.e. harvested) 
power. 

This paper is organised as follows: in Section [TT] the model used for the dynamics of slender structures with 



non- uniform damping in an axial flow is presented. In Section III, the case of uniform damping is considered as a 
reference configuration and the role of fluid loading on destabilisation by damping is discussed. Computations are 
then performed for a neutrally buoyant slender cylinder with moderate fluid loading, and the role of damping on the 



harvested power and flutter response is examined. Section [TV] investigates non-uniform damping distributions and 
seeks optimals on two different families of damping functions, either distributed over the whole structure or focused 
on a particular region. In Section [V] the impact of damping distribution on the flutter dynamics is investigated further 
to understand the fundamental difference between optimal configurations at low and intermediate damping. 



II. FLUID-SOLID MODEL 



We consider a cantilevered (clamped-free, fixed at O) slender structure of length L with crosswise dimension D, 
density p s , stiffness B, and nonuniform structural damping B*(s). The slender solid is immersed in a stream of fluid of 
density p moving at mean speed [Too, and the solid motion is confined to the (e x , e y ) plane (Figure [I]). The equations 
of motion are nondimensionalised by the characteristic system scales: p, L, [Too- 



A. Nonlinear beam model 



The flexible structure is modelled as an inextensible Euler-Bernoulli beam, where r(s) is the position vector in the 
fixed coordinate system (e^e^) and s is the curvilinear coordinate. At each point along the beam, the orientation 
0(s,t) is defined as the angle of the tangent vector t(s) with the horizontal; n(s) is the local normal. The nonlinear 
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equation of motion for the beam subjected to a fluid force, f , is: 
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where the internal tension, z/(s, £), is essentially a Lagrange multiplier to satisfy the inextensibility condition dr/ds = r 
and £(s) = UoqB*(s) / (BL) is the non-dimensional damping distribution. A Kelvin- Voigt damping model is considered 
here, generalising previous contributions [e.g. [TDJ [14j |32j [33] to nonuniform damping distributions. 
The clamped- free boundary conditions must also be satisfied, namely at the fixed end (s = 0): 



and at the free end (s = 1): 
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Consistently with prior work [2TJE5], the nondimensional mass ratio (M*) and flow speed (U*) are defined as: 



M* 
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with A the cross-sectional area of the structure. 



B. Fluid dynamic model 

In the limit of slender structures (D < L), and for purely potential flow upstream of the structure's trailing edge, 
Lighthill's large amplitude elongated-body theory [34] provides a leading order expression for the 'reactive' force P 
applied by the flow on the flapping body, associated with the local transverse motion of each cross-section 

i m a f d(u n rv) d(u n u T n) 1 d(ulr) \ 

M*\ dt ds + 2 ds y K) 

where u T r + u n n = dr/dt — e x is the solid's local velocity relative to the incoming flow, and m a = M a / p s A, where 
M a is the dimensional added mass per unit length of the cross-section. [35] showed that Lighthill's theory compares 
well with Reynolds- averaged Navier-Stokes (RANS) simulations to compute the forces on a swimming fish during 
transient manoeuvres. This reactive force does not include any flow separation associated with the transverse motion 
of each cross-section, and as emphasised in [35], for freely flapping bodies, an additional contribution for the fluid 
force must be included to account for such dissipative effects. Here, an empirical 'resistive' force model is therefore 
added following [36] and [33] : 

r = -l/2C d u n \u n \n, (6) 

where Cd is the empirical drag coefficient, and = 1 is used in the following for circular cross- sections [see|29] for 
a discussion of the impact of this coefficient on the flapping dynamics]. [37] tuned the drag coefficients to show good 
agreement with direct numerical simulations. 

Thus the fluid force, f , in Eq. is modelled as the sum of the reactive (P) and resistive (P) components. In the 
remainder of the paper we assume a neutrally buoyant circular cylinder, therefore p = p s and A = ttD 2 /4. Unless 
otherwise stated, we assume D/L = 0.1. 

Note that the fluid force description is purely local here, and does not explicitly account for wake effects. When the 
slender body assumption is not verified, and in particular, in the case of two-dimensional plates, an explicit description 
becomes necessary [see for example [20] [2TJ [38] . 



C. Energy harvester model 

As noted earlier, energy harvesting is represented as a strain-based damping We focus here on the mean 

nondimensional harvested power: 

V= Jfrrs = ^—2 / t(s)(K 2 ) ds, (7) 

pDLU^ M*U* 2 Jo 
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FIG. 2: (a) Maximum deflection y m &* (solid) and orientation # m ax (dashed) at the free end as a function of the non-dimensional 
flow velocity U* for M* = 12.7. (b) Snapshots of the beam response for U* = 10, 13 and 19 (from top to bottom). 



where & is the mean dimensional harvested power, k is the time derivative of the local curvature k and (•) is the time 
average taken over a period T of the limit-cycle oscillation. V can also be understood as the efficiency of the system. 

As discussed in Section [IJ the presence of structural damping can affect the flutter response. The intensity and 
distribution of damping are characterised by 

£o = / £(s)ds, and i(s) = (8) 
Jo 

Equation ([8| allows us to independently evaluate the impact on the system response of (i) the amount of damping £o 
and (ii) its spatial distribution. 



D. Numerical solution 



Equation ([I]) is solved numerically together with boundary conditions ([2])-([3]) using an iterative second-order implicit 
time-stepping scheme [39 j5 and spatial derivatives are computed using Chebyshev collocation [40]. Conservation of 
energy is ensured by verifying that E = Wf — Q, where 




are respectively the mechanical energy of the system, the dissipated power and the rate of work of the fluid forces, 
and is classically obtained by projecting the equation of motion ([!]) on the solid velocity r and by integrating over 
the entire beam. Also note from (|7|), V = (Q). In the numerical implementation the beam and the flow are initially 
at rest; the flow speed is ramped up to its steady state value and a small perturbation is applied to the vertical flow. 



E. Nonlinear response of an undamped beam 



Prior to analysing the energy harvesting properties of the system, we first examine the undamped flutter response 
of the structure. In Figure [2j we plot the system response for increasing nondimensional flow speed, for a circular 
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FIG. 3: Variation of the critical nutter speed, U c f with damping £o for M* = 2.5 (dotted), 5.1 (dash-dotted), 8.5 (dashed), 12.7 
(solid) and 20 (thick solid) obtained from linear stability analysis for uniform damping. Results of nonlinear computations for 
M* = 12.7 are also presented (squares). 

cylinder (M* ~ 12.7). The critical flow speed at which flutter ensues is verified from linear stability analysis and 
is confirmed with [26 j. Consistent with flutter in plates [33], we note that this instability is a supercritical Hopf 
bifurcation with flow speed. The jumps in the bifurcation curve correspond to the mode switching reported in [41], 
this may also be discerned from the snapshots of the beam at three different flow speeds. For the energy harvesting 
computations performed in the rest of the paper we set the flow speed at U* = 13, corresponding to well-developed 
oscillations and moderate deflections for the undamped configuration. 

III. UNIFORM DAMPING 

Here, a uniform damping distribution £(s) = £o is considered. The dependence of critical flutter speed on the 
damping intensity £o is first analysed for varying values of fluid loading M*. Based on these results, a moderate fluid 
loading is selected to study the power harvesting capacity of the configuration. 

A. Destabilisation by damping: impact on critical flutter speed 

As noted in Section [TJ damping can destabilise flexible structures at sufficiently high fluid loading. [32] shows that 
destabilisation of long pipes with a high mass ratio is associated with a drop in the critical flutter speed. Based on 
this work, we analyse the linear stability of the system and in Figure [3] compare the variation of the critical flutter 
speed, U c f, with damping, £o> a t different values of fluid loading, M*. For lightly loaded structures (M* < 5.5) U c f 
monotonically increases with damping; at higher values (M* > 8.5), U c f decreases with damping until a specific value 
(£ m ) above which U c f increases rapidly; note that £ m increases with M*. 

For the remainder of the energy harvesting analysis, we settle on a value of M* = 12.7 (corresponding to D/L = 
0.1): this choice allows us to analyse the scope for harvesting energy from a system at moderate fluid loading with 
destabilisation by damping. 

B. Harvesting power with constant damping 

Figure [4^ a) presents the evolution of harvested power, V, with damping for 10 -3 < £o < 10. Equation Q becomes 

P= — || K ||i, (11) 

M*TT* 2 
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FIG. 4: (a) Evolution of the harvested power V (solid) and curvature norm || JC ||i (dashed) with the uniform damping intensity 
£o- (b) Distribution of curvature-change JC(s) along the beam for £o = 0, 0.022, 0.22 and 0.47 (from left to right and top to 
bottom, respectively). (M* = 12.7, U* = 13). 

where JC(s) = (k 2 ) and || JC ||i= J* JCds is the Li-norm of tC(s). The evolution of the flutter response with damping 
is examined by plotting the rescaled curvature term, JC = tC/ || /Co ||oo (where JCo(s) is the value of JC for £o = 0) for 
increasing £o (Figure |4)d). At small £o> the response is virtually no different from the undamped case, but for £o > 0.1 
a perceptible change is observed. First, the curvature is redistributed along the entire beam, and as the damping is 
increased further, the response of the beam is damped out globally. This can also be seen from the variation of || JC ||i 
with £o m Figure [4^ a). Notably from ( 

As a result, for small damping, the change in system response is virtually imperceptible from the undamped case, 
and power simply scales linearly with £ - However, damping has a strong effect on the flutter response at larger it 
reduces the amplitude of curvature-change significantly and causes a sharp reduction in V. The strategy to optimally 
harvest power is to find the upper bound on £o below which || JC ||i can be maintained close to (or ideally enhanced 
above) its undamped value. 

IV. NONUNIFORM DAMPING 

The results for constant damping suggest that as long as || JC ||i is maintained at undamped levels, the harvested 
power increases linearly with £o- Figure [4jb) clearly shows that JC(s) varies significantly along the length of the beam. 
Concentrating harvesters around the zone of peak curvature could therefore enhance the harvested power. This idea 
is tested in this section on a reduced functional space for the damping distribution Two families of damping 

functions are considered, namely (a) dispersed and (b) focused damping distributions. Optimisation of the harvested 
power is performed within each family, in anticipation of insights into the global optimal. 

A. Dispersed harvester distribution 

In this section, we are interested in simple non-homogeneous distributions of damping of the form: 

C(s) = Co(l+Ci(s-l/2)) (12) 

characterised by the total damping, 10 -3 < £o < 10, and slope, — 2 < £i < 2. This function family corresponds to 
a dispersed distribution, where the damping is significant over the entire length of the structure. Figure |5|a) shows 
the variation of re scal ed power, V = V/V™* with (£o, where 7^™ ax is the maximum harvested power for constant 
damping (Section |III| ). Figure [5] shows that for all £o> the optimal distribution corresponds to £i =2, when damping 
is distributed increasingly from fixed to free end, and that this linear distribution of damping leads to an increase 



11) it is || JC ||i that directly impacts the harvested power. 
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FIG. 5: Linear distribution (|12|): 

^P/^ ax and (c 



(a) rescaled power contours p/p™ ax for varying £o and £i. 



/C ||i as a function of the damping intensity £o for linearly decreasing (£i = 
dotted) and linearly increasing (£i = 2, solid) damping distributions. (M* = 12.7, C7* = 13). 



(b) Rescaled harvested power 
—2, dashed), constant (£i = 0, 



of the maximum harvested power by 50% compared to the uniform distribution (£i = 0). Through non-uniform but 
simple distributions of damping on the structure, it is indeed possible to enhance the flutter response above that of 
the undamped configuration (Figure^). 



B. Focused harvester distribution 



The results from the reduced order analysis [29] suggest that the optimal distribution ought to be localised at 
specific points on the beam. The chief drawback of the two parameter functions examined in the previous section is 
the inability to consider localised distributions of damping in specific regions. To consider such peaked distributions, 
we now turn to the following three-parameter family of Gaussian damping distribution: 

£00 = £o^Y' t g = e-« a —r, (13) 
II €g 111 

where £o is the total damping in the system, s Q is the centre of the distribution and 1/a is a measure of its spread in 
s. For increasing a energy harvesters are increasingly concentrated around s Q . 

Figure [6] presents the rescaled power V/V™ 3 '* in the (a, s G )-space for small (£o = 0.01) and moderate damping 
(£o = 0.47). For small damping, one sees that power is optimally harvested for dampers focused at s Q = 0.8, the 
position of the maximum of /Co along the beam (Figure ^p). This is quite distinct from the high damping optimal, 
which corresponds to a dispersed distribution (a = 2.7, s Q = 1). Of particular note is that for all £o> a uniform 
distribution (a = 0) is preferred over a concentration of harvesters at s Q < 0.5. 

We next examine a wider range of damping, 10 -3 < £o < 50, and plot the optimal values of (a, s Q ), in Figure 7(a) 



For moderate to large damping, £o > 0.3, a dispersed distribution with peak at s Q = 1 is optimal. Conversely for 
£o < 0.02, harvesters concentrated at s Q w 0.8 (position of peak /C) is the optimal configuration. Note that for these 
computations we set < a < 500, and at these values of £o the upper bound is reach ed; more detailed insights and 
computations on the small damping regime are presented in Section |V| Figure |7(b)| shows the evolution with £ of 
the optimal normalised distribution and illustrates the transition from focused to dispersed profiles when the total 
damping is increased. 
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FIG. 6: Gaussian distribution (13): Maps of rescaled power V /Vj? 1 
(M* 12.7, U* 13). 



for varying a and s Q at (a) £o = 0.01 and (b) £o = 0.47. 
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FIG. 7: Gaussian function optimal configuration: (a) Spread parameter, a, (open circles) and centre location, s G , (filled 
boxes) corresponding to the optimal Gaussian distribution for a given £o- (b) Rescaled damping distribution at discrete 
£o = 0.004, 0.025, 0.25, 0.47 (solid, dashed, dashed-dot and dotted curves, respectively, indicated on (a) with vertical lines of 
corresponding description); superimposed is the linear optimal distribution (£i = 2; thick solid curve). (M* = 12.7, U* = 13). 



It is possible to determine the optimal harvested power for each value of £o; in Figure [8] these optimals are compared 
with the results for uniform and linear distributions. Strikingly, the peak power values are coincident for linear 
and Gaussian distributions, which is consistent with the similarity in distributions (Figure 7(b)). Because of the 
fundamentally different structures of the distribution used, this result suggests that the optimal obtained with such 
simple functions represents a good approximation of the absolute optimal. Departing from this maximum of harvested 
power, in particular at small damping we see that concentrated harvesting is superior to the optimal linear distribution. 

These results confirm that nonuniform damping distributions can be advantageously employed to enhance the power 
harvesting capacity. Focused damping distributions are optimal for small damping, whereas dispersed distributions 
are preferential at higher damping. The following section investigates in more details this transition by considering 
the impact of damping on the nonlinear dynamical response of the system. 
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FIG. 8: Optimal power values corresponding to Gaussian optimal configuration in Figure [7] (squares) compared to the optimal 
power obtained for linear distribution (£i = 2; thick curve) and constant distribution (£i = 0; dotted curve). 



V. DISCUSSION: IMPACT OF LOCALISED DAMPING 



The above computations suggest a very different impact of damping on the dynamics of the structure depending 
on the magnitude of £o? thereby leading to quite different optimal strategies to maximise the harvested power. In this 
section, we first consider the limit of asymptotically small damping before studying the effect of focused dampers on 
the body's dynamics. 



A. Optimal distribution for asymptotically small damping 

For asymptotically small damping (|| £ Hoo^ 1, with || £ the maximum value of £(s) for s G [0 1]), the flapping 
dynamics is not modified at leading order so that /C = /Co + O(||£|| 00 ). Thus 

v ~ Wu^ jf miCo ds ~ 11 K ° l|co (14) 



and this upper bound can be approached asymptotically using, for example, the Gaussian distribution (13) centred 
on the maximum of /Co and increasing a. More precisely, any peaked distribution of damping of Li-norm £o centred 
on this maximum and with a typical width (here a -1 / 2 ) much smaller than the length scale A associated with the /Co 
peak width, should approach the theoretical upper bound V th = £o || ^oljpo /(M*Z7* 2 ), provided that the maximum 
damping remains small enough that the approximation for V in Eq. (fiS still holds (which in the present case is 
equivalent to keeping a 1 / 2 ^ bounded). 

This conjecture is tested numerically using a Gaussian distribution, Eq. ( [13] ^ with s Q = 0.8, and by computing the 
harvested power for small damping: 2.2 x 10 -5 < £o < 2.2 x 10 -2 . Figure"^ shows P/P th and confirms that this 
conjecture indeed holds since for increasing a, the power asymptotically approaches its optimal for small enough £ - 
When £o > 2.2 x 10 -3 , the asymptotic value V th can not be reached any more: as a is increased, the effect of the 
maximum damping a^o on the flutter dynamics is important even before the dampers are focused enough for V to 
approach V th . 

Two essential results are illustrated here. For sufficiently small damping, the nonlinear response of the structure 
remains unchanged so the best strategy is to focus all the harvesters in the region of maximum curvature-change. 
However finite levels of damping significantly modify the nonlinear dynamics of the beam, so at these levels a narrow 
focusing of damping is sub-optimal; next we investigate further the behaviour at finite damping. 
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FIG. 9: Focused harvesters at s Q — 
harvesters for small damping values: 
12.7, U* = 13). 



0.81 for small damping: Curves show the normalised power for increasingly focused 
= 2.2 x 10 -5 : 10 -2 (solid, dashed, dash-dot, dotted curves, respectively). (M* = 



B. Impact of finite and localised damping 



The computations from Section [TV] show that for Gaussian distributions and £o > 0.02 a defocusing of harvesters is 
preferable, as the system response becomes increasingly influenced by damping. In order to understand the impact on 
the system's response, the Gaussian distribution (13) is used in the high damping range and a is varied over a range 
that transitions the distribution from dispersed to focused (0 < a < 10 2 ) with £o = 0.47 and s Q = 0.45 (corresponding 
to the positio n of || K ||oo for a = 0). Consistent with Figure [6j the harvested power in Figure [lb (a) | decreases with 
a, and Figure [To(b)| shows the evolution of JC(s). 

For dispersed distributions (a < 10), JC(s Q ) decreases with a and the zone of maximum curvature-change shifts to 
regions of reduced damping. Nonetheless nonzero damping in this zone is adequate to retrieve some of the energy. 
For a > 10 damping becomes increasingly focused and the flexible body does not deform anymore near the damper 
(JC(sq) ~ 0): no energy is harvested anymore since no deformation occurs near the harvester's position, and this is 
reflected in the sharp drop in power in Figure |10(b)| This illustrates the main effect of a focused distribution for 
finite damping: we see a redistribution of the solid's deformation to regions with little or no damping. This is also 
the reason why a focused damping distribution is not appropriate for optimal energy harvesting in the finite damping 
range. 



VI. CONCLUSIONS 



In this work, we considered the possibility of harvesting energy from a slender body fluttering in an axial flow, and 
in particular the impact of harvester-distribution on the performance of the system as well as potential optimisation 
strategies. To this end, a simplified fluid-solid model was proposed with energy harvesting represented as nonuniform 
structural damping. Depending on the fluid-to-solid inertia ratio, damping can actually enhance the flutter response 
of the structure, as well as reduce the critical flow velocity above which the system can operate. 

For uniform damping distributions, maximising the harvested power appears as a trade-off on the damping intensity: 
for small damping, the flutter response is only weakly modified and the harvested power increases as more damping 
is added to the system. For higher damping however, the dynamical response can be strongly modified and the self- 
sustained oscillations are eventually mitigated. The effect of a nonuniform distribution of harvesters along the structure 
was considered next. We show that even simple nonuniform distributions such as linear and Gaussian functions, can 
lead to an increase in harvested peak power on the order of 50%. The similarity in the optimal distribution and 
performance obtained through an optimisation on two fundamentally different families of distributions suggests that 
simple strategies can capture rather well the characteristics of the global optimal distribution. 

Investigating further into the relationship between damping and flutter response, we showed that for small damping, 
localised harvesting is optimal as it takes full advantage of the system's maximum curvature without impacting its 
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FIG. 10: Focused harvesters at s Q = 0.45 for large damping (£ = 0.47): (a) power dependence on a and (b) JC(s) for 
a = 0, 1, 10, 10 2 (solid, dashed, dashed-dot, dotted curves, respectively). (M* = 12.7, U* = 13). 



dynamics significantly. On the other hand, for finite damping, focused distributions perform rather poorly as the 
beam response adapts to rigidify the damped region, leading to negligible harvested power. Instead, nonuniform 
distributed damping over the entire length of the system becomes optimal. 

Returning to the result from the simple bi-articulated model [29] , we note a clear difference in the optimal con- 
figurations. Whilst the reduced order model optimal has dampers focused at the fixed end, the continuous optimal 
requires a dispersed distribution with minimal damping at the fixed end and increasing to a maximum at the free 
end. However, we find a crucial difference between the two system configurations: whilst the bi-articulated system 
has a single moving region of curvature (the second articulation) that is responsible for driving the instability, defor- 
mations may occur all along the length of the beam in the continuous system. Therefore in the latter, curvature can 
be displaced away from regions with focused damping while still maintaining the flutter dynamics. In both cases, a 
careful understanding of the dynamics of the system is necessary to determine the optimal nonuniform positioning of 
energy harvesters. 
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